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Abstract 

In the framework of the Constrained Minimal Supersymmetric Standard Model 
(CMSSM) we discuss the impact of the pseudo-scalar Higgs boson in delineating 
regions of the parameters which are consistent with cosmological data and E821 
data on the anomalous magnetic moment of muon. For large values of the param- 
eter tan/3 > 50, cosmologically allowed corridors of large mo, M 1 / 2 are opened, 
due to the s-channel pseudo-scalar exchange in the pair annihilation of the lightest 
of the neutralinos to bb or rf , which dominates in this region. However, no such 
corridors are found for values tan (3 < 50. Combining cosmological and E821 data 
puts severe upper limits on sparticle masses. We find that at LHC, but even at a 
e + e~ linear collider with center of mass energy y/s = 800 GeV, such as TESLA, 
supersymmetry can be discovered, if it is based on the CMSSM. 



1 Introduction 



In the framework of supersymmetric models with P-parity conservation, it has been 
argued that for large tan f3 the neutralino relic density (0,% h%) can be compatible with 
the recent cosmological data which favour small values for fly. h^. In this regime the 
neutralino (x) pair annihilation through s-channel pseudo-scalar Higgs boson (A) ex- 
change leads to an enhanced annihilation cross sections reducing significantly the relic 
density 0, while the heavy CP-even Higgs (H) exchange is P-wave suppressed and 
not that important. The importance of this mechanism, in conjunction with the recent 
cosmological data which favour small values of the Dark Matter (DM) relic density, has 
been stressed in |2[] . The same mechanism has been also invoked in Ref . || and [[| where 
it has been shown that it enlarges the cosmologically allowed regions. In fact Cosmology 
does not put severe upper bounds on sparticle masses, and soft masses can be in the TeV 
region, pushing up the sparticle mass spectrum to regions that might escape detection 
in future planned accelerators. Such upper bounds are imposed, however, by the recent 
g — 2 E281 data || constraining the CMSSM in such a way that supersymmetry will be 
accessible to LHC or other planned e + e~ linear colliders if their center of mass energy 
is larger than about 1.2 TeV ||. The bounds put by the g — 2 has been the subject of 
intense phenomenological study @-^|, and although the situation has not been definitely 
settled, supersymmetry emerges as a prominent candidate in explaining the discrepancy 
between the Standard Model predictions and experimental measurements. 

In this study we undertake the problem of calculating the neutralino relic density, 
in the framework of the CMSSM, paying special attention to the pseudo-scalar Higgs 
exchange which dominates in the large tan/3 region. In this regime the cosmologically 
allowed domains depend sensitively on this mechanism and in conjunction with the 
bounds put by the g — 2 measurements can severely constrain the CMSSM predictions. 
In particular, we find that cosmologically allowed corridors of large mo, M1/2 values open 
up for tan (3 > 50, which however have little overlap with the regions allowed by the E821 
data. The constraints imposed on the sparticle spectrum and the potential of discovering 
CMSSM in future accelerators are discussed. 
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2 The role of the pseudo-scalar Higgs boson mass 

The XX fusion to the pseudo-scalar Higgs boson, A, which subsequently decays to a bb or 
a tt, becomes the dominant annihilation mechanism for large tan/?, when the pseudo- 
scalar mass rriA approaches twice the neutralino mass, ~ 2m^. In fact by increasing 
tan/? the mass tjia decreases, while the neutralino mass remains almost constant, if 
the other parameters are kept fixed. Thus vha is expected eventually to enter into the 
regime in which it is close to the pole value rriA = 2m^, and the pseudo-scalar Higgs 
exchange dominates. It is interesting to point out that in a previous analysis of the direct 



DM searches [10], we had stressed that the contribution of the CP-even Higgs bosons 
exchange to the LSP-nucleon scattering cross sections increases with tan (3. Therefore in 
the large tan (3 regime one obtains the highest possible rates for the direct DM searches. 
Similar results are presented in Ref. |P"T|| . 

In the framework of the CMSSM the chargino mass bound as well as the recent LEP 



Higgs mass bound [ff2| already exclude regions in which \ h as a large Higgsino compo- 
nent, and thus in the regions of interest the x is mainly a bino. A bino is characterized 
by a very small coupling to the pseudo-scalar Higgs A, however the largeness of tan (3 
balances the smallness of its coupling giving a sizeable effect when ttla — 2m^, making 
the s-channel pseudo-scalar exchange important. 

It becomes obvious from the previous discussion that an unambiguous and reliable 
determination of the v4-mass, itla, is demanded before one embarks on to calculate the 
neutralino relic density especially in regions where the s-channel pseudo-scalar exchange 
dominates. In the constrained SUSY models, such as the CMSSM, rriA is not a free 
parameter but is determined once m , M1/2, A as well as tan/? and the sign of /x, 
sign(fi), are given. rriA depends sensitively on the Higgs mixing parameter, mf, which 
is determined from minimizing the one-loop corrected effective potential. A subtlety 
arises for large tan/3 values since the corrections are relatively large due mainly to the 
smallness of the Higgs mixing parameter. In order to handle this we calculate the 



m,-, m 



ti""t 2 



effective potential using as reference scale the average stop scale Qi ~ ^ 
At this scale the contributions of the third generation sfermions are small. However, 
other contributions may not be negligible at this scale and should be properly taken 
into account. In particular, the neutralino and chargino contributions to the effective 
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potential should be included for a more reliable calculation. These do not vanish at Qi 
since their masses, determined by the gaugino masses M\,Mi and the \x value, may be 
quite different from Q~ t . Their inclusion to the effective potential improves the mass of 
the A-Higgs, as this is calculated from the effective potential, yielding a result that is 
scale independent and approximates the pole mass to better than 2% if we also include 
the scale dependent logarithmic contributions from the wave function renormalization 
n(0) — U(m\) |f[4j| . I the present work and for a more accurate determination of the 



pseudo-scalar Higgs we use the pole mass using the expression of Ref . JT5[ . 

A more significant correction, which drastically affects the pseudo-scalar mass arises 
from the gluino-sbottom and chargino-stop corrections to the bottom quark Yukawa 



coupling ||15|-|20|j. Taking these effects into account the tree-level relation between the 
bottom mass and the corresponding Yukawa coupling is modified according to 



nib = vi(hb + Ahb tan (3 



(1) 



where Ahb is 
Ah b 
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(2) 



In this the first (second) term is the gluino, sbottom (chargino, stop) corrections. In the 
second term we have ignored the small electroweak mixing effects. The function G in 
Eq. (D is the one used in Refs. ]l7|,f[^,^0| which in [JL7]] is denoted by I. It is known that 
the proper resummation of these corrections is important for a correct determination of 
hb [ l9|,|20| . Eqs. ([!]) and (@) agree with those of Ref. [20| and therefore these corrections 
have been properly resummed []. 

These are very important, especially the SQCD corrections, and should be duly taken 
into account. The important point is that these affect differently the /i > and fi < 
values of itla in a drastic way. In figure [I] we show the behaviour for the pseudo-scalar 
Higgs boson mass as function of the parameter tan f3 for the inputs shown in the figure, 
for both signs of fi, with and without the aforementioned corrections. Also plotted is 
the value of the quantity 2 m^. One observes that in the absence of these correction to 
the bottom Yukawa coupling, the pseudo-scalar Higgs mass in the two cases differ little 



1 Our fi and the soft gaugino masses differ in sign from those of Ref. 117,201, while At has the same 



sign. 
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and meet the 2 line at tan/3 ~ 40. Obviously the tan/? region around this value 
leads to enhanced neutralino annihilation cross sections through ^4-exchange since we 
are in the vicinity of a pole. However, when the aforementioned corrections are taken 
into account the values of the pseudo-scalar mass in the two cases split as shown in 
the figure. That corresponding to the /i > is moving upwards and that to /i < 
downwards. Thus only the second can reach the pole value 2 at a smaller tan (3 
however. The mass corresponding to > stays away from this, never reaching it, 
at least in the case shown, since for higher values of tan/3 we enter regions which are 
theoretically forbidden. Actually in these regions Electroweak Symmetry Breaking does 
not occur. The fi < case does not stay comfortably well with the b — > s + 7 process, 
as well as with the observed discrepancy of the g — 2 data, if the latter are attributed to 
supersymmetry, and therefore we shall discard it in the sequel. 

It is obvious from the previous discussion that the crucial parameter in this analysis 
is the ratio m^/^m^. The calculation of the pseudo-scalar Higgs boson mass rriA we 
have discussed in detail before. For the calculation of the x mass we use the one-loop 
corrections of Ref. [I5|. These result to corrections as large as 5%, in some cases, reducing 
the aforementioned ratio, making it to be closer to the pole value. 



3 Numerical results— Discussion 

Before embarking on to present our results we shall comment on our numerical analysis 
employed in this paper. This will be useful when comparing the results of this paper 
with those of other works, which use different numerical schemes in determining the 
mass parameters of the CMSSM. The predictions for the sparticle spectrum, including 
the mass of the pseudo-scalar Higgs, as well as the calculation of the relic density itself, 
may be sensitive to some of the parameters and the particular scheme employed. 

In our analysis we use two- loop renormalization group equations (RGE), in the DR 
scheme, for all masses and couplings involved. The unification scale Mqut is defined as 
the point at which the gauge couplings &i and 0.2 meet but we do not enforce unification 
of the strong coupling constant 0J3 with 0:1,2 at Mqut- The experimental value of the 
MS strong coupling constant at Mz is an input in our scheme and this is related to 
&3 through a s {Mz) = &3(Mz)/(l — Ads), where Ad 3 are the threshold corrections. 
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Enforcing unification of A3 with the rest of the gauge couplings usually results in values 
for a s (Mz) larger than the experimental values, if the two-loop RGE's are used. For 
this reason we abandon gauge coupling unification. For the determination of the gauge 
couplings ai t 2 we use as inputs the electromagnetic coupling constant qj , the value of 
the Fermi coupling constant Gp, and the Z-boson mass Mz- From these we determine 
the weak mixing angle through the relation s c = na I V2M 2 Z G F {1 - Af). The value 



of the electromagnetic coupling constant at Mz in the DR scheme is calculated through 
&(Mz) = Qq/( 1 — Ad em ), where Aa em are the appropriate threshold corrections (for 
details see Ref. ||15||). In each iteration s 2 and a are extracted and from these the values 
of ai t 2 at Mz are determined. In the equations above all hatted quantities are meant to 



be in the DR scheme. 

Our remaining inputs, in running the RGE's, are as usual the soft SUSY breaking 
parameters m , M1/2, Aq, tan (3 and the sign of the parameter \x. The top and tau physical 



masses, M t ,M T , as well as the MS bottom running mass mb(rrib) are also inputs. The 
default values for the aforementioned masses are M t = 175 GeV, M T = 1.777 GeV and 
rribimb) = 4.25 GeV. For the determination of the bottom and tau running masses, and 
hence their corresponding Yukawa couplings, at the scale Mz, we run SU C (3) x U em (l) 
RGE's using three-loop RGE's for the strong coupling constant. We also include two- 
loop contributions in the electromagnetic coupling and mixings of the electromagnetic 
with the strong coupling constant. The latter is as important as the three- loop strong 



coupling constant contribution to the RGE's. At the end of the running MS masses 



are converted to DR in the usual way. This determines the bottom and tau Yukawa 
couplings at Mz- We recall that the corrections to the bottom Yukawa coupling of 
Eq. (0) should be duly taken into account. For the determination of the top Yukawa 
coupling at M t we take into account all dominant corrections relating the pole to its 
running mass. By running the RGE's we can have the top Yukawa coupling at Mz- 



Thus our analysis resembles that followed in Ref. flT5 



For the determination of the Higgs and Higgsino mixing parameters, m\ and /1, we 
solve the minimization conditions with the one-loop corrected effective potential in which 
all particle contributions are taken into account. The minimization is performed using 
as reference scale the average stop scale Qi ~ y/ m ti m t 2 • Thus in each run we determine 
ml(Qi),^(Qt). 
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Regarding the calculation of the lightest supersymmetric particle (LSP) relic abun- 
dance, the Boltzmann equation is solved numerically using the machinery outlined in 
Ref. [0. The coannihilation effects in regions where the right-handed stau, tr, ap- 
proaches in mass the LSP, are properly taken into account. Caution should be taken in 
regions where the ratio rriA/2m^ is close to unity. This signals the vicinity of a pole in 
which case the traditional non-relativistic expansion breaks down. On the pole the an- 
nihilation cross section through the pseudo-scalar Higgs s-channel exchange is large and 
its width is important in determining its size. The rescaled pseudo-scalar Higgs boson 
width is r a/ 111 a — 10~ 2 , that is resembles that of the Z-boson, and hence for values of 
the ratio m^/^m^ larger than about 1.2 (see pl| ) we are away from the pole region. 
We have found that in the mo,Mi/ 2 plane, and for both parameters less than 1 TeV, 
this ratio approaches unity only for very large values of tan/3 > 50. For such values of 
tan /3 the pseudo-scalar Higgs dominates the x P arr annihilation, leading to large cross 
sections and therefore cosmologically acceptable relic densities. Thus cosmologically al- 
lowed mo, Mx/2 corridors open up for tan/3 > 50 which were absent for lower values of 
tan p. These are the same corridors observed in the analysis of Ref. || which however 
show up for lower values of tan p. 

In the panels shown in figure we display our results by drawing the cosmologically 

allowed region 0.08 < fl^ < 0.18 (dark green) in the m ,Mi/2 plane for values of 

tan/3 equal to 40, 45, 50 and 55 respectively. Also drawn (light green) is the region 

0.18 < h,Q < 0.30. In the figures shown the default values for the top, tau and bottom 

masses are assumed. The remaining inputs are shown on the top of each panel. The 

solid red mark the region within which the supersymmetric contribution to the anomalous 

magnetic moment of the muon falls within the E821 range a^ USY = (43.0 ±16.0) x 10~ 10 . 

The dashed red line marks the boundary of the region when the more relaxed lower bound 

11.2 x 10 -10 < a^ USY is used [[J, corresponding to the 2cr lower bound of the E821 range. 

Along the blue dashed-dotted lines the light CP-even Higgs mass takes values 113.5 GeV 

(left) and 117.0 GeV (right) respectively. The line on the left marks therefore the recent 

LEP bound on the Higgs mass (T^|. Also shown is the chargino mass bound 104 GeV Q. 

The shaded area (in red) at the bottom of each figure, labelled by TH, is theoretically 

2 In the context of our analysis focus point regions || show up for smaller values of the top mass. 
At this point we therefore agree with the findings of Ref. In any case the bulk of the focus point 

region appears for rather large values of toq and hence they are not favoured by the g — 2 data. 
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disallowed since the light stau is lighter than the lightest of the neutralinos. From the 
displayed figures we observe that for values of tan (3 up to 50 the cosmological data put an 
upper bound on the parameter m®. However, there is practically no such upper bound 
for the parameter My 2, due to the coannihilation effects M which allow for Mi/ 2 as 
large as 1700 GeV within the narrow coannihilation band lying above the theoretically 
disallowed region. For tan (3 = 55 a large region opens up within which the relic density 
is cosmologically allowed. This is due to the pair annihilation of the neutralinos through 
the pseudo-scalar Higgs exchange in the s-channel. As explained before, for such high 
tan/3 the ratio m^/2m^ approaches unity and the pseudo-scalar exchange dominates 
yielding large cross sections and hence small neutralino relic densities. In this case the 
lower bound put by the g — 2 data cuts the cosmologically allowed region which would 
otherwise allow for very large values of mo,Mi/ 2 . The importance of these corridors 
has been stressed in the analysis of M and ||. However, in our case these show up at 
much higher values of the parameter tan (3. We should remark at this point that in our 
analysis we use the value of a s t rong (Mz) as input and relax unification of the a 3 gauge 
coupling with the others. For reasons already explained, in the constrained scenario it 
is almost impossible to reconcile gauge coupling unification with a value for a s t ro „ 9 (M z ) 
consistent with experiment due to the low energy threshold effects. This change affects 
drastically the values of other parameters and especially that of the Higgsino (/i) and 
Higgs mixing parameters that in turn affect the pseudo-scalar Higgs boson mass 
which plays a dominant role. For the tan f3 = 55 case, close the highest possible value, and 
considering the conservative lower bound on the muon's anomalous magnetic moment 
a susY > -q 2 x 10~ 10 and values of Q% h% in the range 0.13 ±0.05, we find that the point 
with the highest value of m is (in GeV) at (mo,M 1 / 2 ) = (950,300) and that with the 
highest value of My 2 is at (m , M1/2) = (600, 750). The latter marks the lower end of the 
line segment of the boundary a^ USY = 11.2 x 10~ 10 which amputates the cosmologically 
allowed stripe. For the case displayed in the bottom right panel of the figure |2| the upper 
mass limits put on the LSP, and the lightest of the charginos, stops and the staus are 
rriy. < 287, < 539, mi < 1161, < 621 (in GeV). Allowing for A 7^ values, the 
upper bounds put on m , Mi/ 2 increase a little and so do the aforementioned bounds 
on the sparticle masses. Thus it appears that the prospects of discovering CMSSM at 
a e + e~ collider with center of mass energy ^Js = 800 GeV, such as TESLA, are not 
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guaranteed. However in the allowed regions the next to the lightest neutralino, has 
a mass very close to the lightest of the charginos and hence the process e + e~ — > xx'i 
with x' subsequently decaying to x + or X + 2 jets, is kinematically allowed for such 
large tan/3, provided the energy is increased to at least yfs = 900 GeV. It should be 
noted however that this channel proceeds via the t-channel exchange of a selectron is 
suppressed due to the heaviness of the exchanged sfermion. 

The situation changes, however, when the strict E821 limits are imposed a^ USY = 
(43.0 ± 16.0) x 10~ 10 . For instance in the tan (3 = 55 case displayed in figure ^| there is no 
cosmologically allowed region which obeys this bound. For the other cases, tan/3 < 50, 
the maximum allowed M1/2 is about 475 GeV, occurring at mo — 375 GeV, and the 
maximum mo is 600 GeV when M1/2 — 300 GeV. The upper limits on the masses of 
the sparticles quoted previously reduce to < 192, < 353, m$ < 775, m f < 436 
all in GeV. However, these values refer to the limiting case A = 0. Scanning the 
parameter space allowing also for A 7^ we obtain the upper limits displayed in the 
table [I]. In this the unbracketed values correspond to the E821 limits on the g — 2. For 
completeness we also display, within brackets, the bounds obtained when the weaker 
lower bound a^ USY > 11.2 x 10~ 10 is imposed. We see that even at TESLA with center 
of mass energy y/s = 800 GeV, the prospects of discovering CMSSM are guaranteed in 
the e + e~ — > x + XT if the E821 bounds are imposed. 

In the figure ||| we display in the (M1/2, mo) plane the points which are consistent both 
with the muon's anomalous magnetic moment bounds mentioned before and cosmology, 
as well as with the other accelerators data. Each of the points is taken from a sample 
of 40,000 random points in the part of the parameter space defined by m < 1.5 TeV, 
M1/2 < 1-5 TeV, \A \ < 1 TeV and 2 < tan/3 < 55. All the points are consistent with 
the cosmological bound £1% = 0.13 ± 0.05. The plus points (colored in blue) are those 
consistent with the E821 bound 27 x 10~ 10 < af, USY < 59 x 10~ 10 , while the diamonds 
(colored in green) are consistent with the more relaxed bound 11.2 x 10~ 10 < a^ USY < 
59 x 10 -10 . The points are grouped in regions, separated by dashed contours, each of 
which constitutes the boundary of tan/3 with the value shown beneath. In the region 
designated as tan/3 = 55 all points have 55 > tan/3 > 50. It is seen clearly that only 
a few points in the tan/3 > 50 case can survive the E821 bound. For tan/3 < 50 the 
parameter M1/2 cannot be larger than about 500 GeV, attaining its maximum value at 
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mo — 400 GeV , and the maximum mo is 725 GeV occurring at M\/2 — 275 GeV. The 
upper limits put on mo, My 2 result to the sparticle mass bounds displayed in the table |I[ 

4 Conclusions 

In this work we have undertaken the problem of calculating the neutralino relic density 
in the framework of the CMSSM, by paying special attention to the pseudo-scalar Higgs 
exchange mechanism which is dominant in the large tan f3 region. Imposing the bounds 
a u U = (43.0 ± 16.0) x 10 -10 on the muon's anomalous magnetic moment, put by the 
BNL E821 experiment, in combination with the cosmological data Q% = 0.13 ± 0.05, 
severely restricts the sparticle spectrum. We found that the pseudo-scalar Higgs exchange 
mechanism opens cosmo logically allowed corridors, of high mo, M\/2, only for very large 
tan j3 > 50, which, however, have little overlap with the regions allowed by the E821 data. 
In fact only a few isolated points in the parameter space with tan (3 > 50 can survive 
the restrictions imposed by both data. The bounds put on the sparticle spectrum can 
guarantee that in LHC but also in a e + e~ linear collider with center of mass energy 
yfs = 800 GeV, such as TESLA, CMSSM can be discovered. The guarantee for a 
TESLA machine with this energy is lost in a charged sparticle final state channel, if the 
lower bound on the value of g — 2 is lowered to its ~ 2a value, but not for the LHC. In 
this case only by increasing the center of mass energy to be ~ 1.2 TeV, a e + e~ linear 
collider can find CMSSM in f f * or x + X~ channels. 
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184 (306) 
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376 (686) 
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20 


154 (255) 


268 (457) 


175 (274) 


603 (990) 


116 (118) 


30 


191 (310) 


338 (560) 


212 (312) 


740 (1200) 


117 (118) 


40 


201 (340) 


357 (617) 


274 (353) 


785 (1314) 


117 (119) 


50 


208 (357) 


371 (646) 


440 (427) 


822 (1357) 


117 (119) 


55 


146 (311) 


260 (563) 


424 (676) 


606 (1237) 


115 (117) 



Table 1: Upper bounds, in GeV, on the masses of the lightest of the neutralinos, 
charginos, staus, stops and Higgs bosons for various values of tan/? if the the E821 
bounds are imposed. The values within brackets represent the same situation when the 
weaker bounds 11.2 x 10~ 10 < a SUSY < 59.0 x 10~ 10 are used (see main text). 
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Figure 1: The pseudo-scalar Higgs masses for \i > and fi < as function of tan/? (solid 
lines). The dashed lines are the same masses when the supersymmetric corrections to 
the bottom Yukawa coupling are ignored. 
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A =0, tan/3=55, fx>0 (m = 4.25 GeV) 
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Figure 2: Cosmo logically allowed regions of the relic density for four different values of 
tan p in the (M1/2, mo) plane. The remaining inputs are shown in each figure. The mass 
of the top is taken 175 GeV. In the dark green shaded area 0.08 < < 0.18. In the 

light green shaded area 0.18 < fl^ < 0.30 . The solid red lines mark the region within 
which the supersymmetric contribution to the anomalous magnetic moment of the muon 



is a 



SUSY 



(43.0 ± 16.0) x 10 . The dashed red line is the boundary of the region for 
which the lower bound is moved to 11.2 < 10 10 af, USY . The dashed-dotted blue lines are 
the boundaries of the region 113.5 GeV < muiggs < 117.0 GeV. 
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/x>0 
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Figure 3: In the (Mi/ 2 ,m ) plane, we display all points compatible with a^ USY = 
(43.0 ± 16.0) x 10~ 10 (+) and 11.2 x 10~ 10 < a s ^ USY < 59.0 x 10~ 10 (o). All the points 
are consistent with the cosmological bound = 0.13 ± 0.05 and they are grouped in 

regions, separated by dashed contours each of which is the boundary of tan (3 with the 
value shown beneath. In the top region, designated by tan/5 = 55, the parameter tan/? 
takes values between 50 and 55. 
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